A novel multi-variate immunological approach, reveals immune variation associated with environmental conditions, and co-infection in the koala (Phascolarctos cinereus).

External signs of disease are frequently used as indicators of disease susceptibility. However, immune profiling can be a more effective indicator to understand how host responses to infection may be shaped by host, pathogen and environmental factors. To better inform wildlife health assessment and research directions, we investigated the utility of a novel multivariate immunophenotyping approach examining innate and adaptive immune responses in differing climatic, pathogen co-infection and demographic contexts across two koala (Phascolarctos cinereus) populations in New South Wales: the Liverpool Plains (LP), and Southern Highlands to South-west Sydney (SHSWS). Relative to the comparatively healthy SHSWS, the LP had greater and more variable innate immune gene expression (IL-1ß, IL-6), and KoRV transcription. During extreme heat and drought, koalas from the LP displayed upregulation of a stress pathway gene and reduced adaptive immune genes expression, haematocrit and plasma protein, suggesting the possibility of environmental impacts through multiple pathways. In those koalas, KoRV transcription status, Chlamydia pecorum infection loads, and visible urogenital inflammation were not associated with immune variation, suggesting that immune markers were more sensitive indicators of real-time impacts than observed disease outcomes.

Contrasting segregation patterns among endogenous retroviruses across the koala population.

Koalas (Phascolarctos cinereus) have experienced a history of retroviral epidemics leaving their trace as heritable endogenous retroviruses (ERVs) in their genomes. A recently identified ERV lineage, named phaCin-ß, shows a pattern of recent, possibly current, activity with high insertional polymorphism in the population. Here, we investigate geographic patterns of three focal ERV lineages of increasing estimated ages, from the koala retrovirus (KoRV) to phaCin-ß and to phaCin-ß-like, using the whole-genome sequencing of 430 koalas from the Koala Genome Survey. Thousands of ERV loci were found across the population, with contrasting patterns of polymorphism. Northern individuals had thousands of KoRV integrations and hundreds of phaCin-ß ERVs. In contrast, southern individuals had higher phaCin-ß frequencies, possibly reflecting more recent activity and a founder effect. Overall, our findings suggest high ERV burden in koalas, reflecting historic retrovirus-host interactions. Importantly, the ERV catalogue supplies improved markers for conservation genetics in this endangered species.

Genetic markers of Chlamydia pecorum virulence in ruminants support short term host-pathogen evolutionary relationships in the koala, Phascolarctos cinereus.

In ruminants infected with Chlamydia pecorum, shorter lengths of coding tandem repeats (CTR) within two genes, the inclusion membrane protein (incA) and Type III secretor protein (ORF663), have been previously associated with pathogenic outcomes. In other chlamydial species, the presence of a chlamydial plasmid has been linked to heightened virulence, and the plasmid is not ubiquitous in C. pecorum across the koala's range. We therefore investigated these three markers: incA, ORF663 and C. pecorum plasmid, as potential indicators of virulence in two koala populations in New South Wales with differing expression of urogenital chlamydiosis; the Liverpool Plains and one across the Southern Highlands and South-west Sydney (SHSWS). We also investigated the diversity of these loci within strains characterised by the national multi-locus sequence typing (MLST) scheme. Although CTR lengths of incA and ORF663 varied across the populations, they occurred only within previously described pathogenic ranges for ruminants. This suggests a relatively short-term host-pathogen co-evolution within koalas and limits the utility of CTR lengths for incA and ORF663 as virulence markers in the species. However, in contrast to reports of evolution of C. pecorum towards lower virulence, as indicated by longer CTR lengths in ruminants and swine, CTR lengths for ORF663 appeared to be diverging towards less common shorter CTR lengths within strains recently introduced to koalas in the Liverpool Plains. We detected the plasmid across 90% and 92% of samples in the Liverpool Plains and SHSWS respectively, limiting its utility as an indicator of virulence. It would be valuable to examine the CTR lengths of these loci across koala populations nationally. Investigation of other hypervariable loci may elucidate the evolutionary trajectory of virulence in C. pecorum induced disease in koalas. Profiling of virulent strains will be important in risk assessments for strain movement to naïve or susceptible populations through translocations and wildlife corridor construction.

Characterisation of the koala (Phascolarctos cinereus) pouch microbiota in a captive population reveals a dysbiotic compositional profile associated with neonatal mortality.

BACKGROUND: Captive koala breeding programmes are essential for long-term species management. However, breeding efficacy is frequently impacted by high neonatal mortality rates in otherwise healthy females. Loss of pouch young typically occurs during early lactation without prior complications during parturition and is often attributed to bacterial infection. While these infections are thought to originate from the maternal pouch, little is known about the microbial composition of koala pouches. As such, we characterised the koala pouch microbiome across the reproductive cycle and identified bacteria associated with mortality in a cohort of 39 captive animals housed at two facilities. RESULTS: Using 16S rRNA gene amplicon sequencing, we observed significant changes in pouch bacterial composition and diversity between reproductive time points, with the lowest diversity observed following parturition (Shannon entropy - 2.46). Of the 39 koalas initially sampled, 17 were successfully bred, after which seven animals lost pouch young (overall mortality rate - 41.18%). Compared to successful breeder pouches, which were largely dominated by Muribaculaceae (phylum - Bacteroidetes), unsuccessful breeder pouches exhibited persistent Enterobacteriaceae (phylum - Proteobacteria) dominance from early lactation until mortality occurred. We identified two species, Pluralibacter gergoviae and Klebsiella pneumoniae, which were associated with poor reproductive outcomes. In vitro antibiotic susceptibility testing identified resistance in both isolates to several antibiotics commonly used in koalas, with the former being multidrug resistant. CONCLUSIONS: This study represents the first cultivation-independent characterisation of the koala pouch microbiota, and the first such investigation in marsupials associated with reproductive outcomes. Overall, our findings provide evidence that overgrowth of pathogenic organisms in the pouch during early development is associated with neonatal mortality in captive koalas. Our identification of previously unreported, multidrug resistant P. gergoviae strains linked to mortality also underscores the need for improved screening and monitoring procedures aimed at minimising neonatal mortality in future. Video Abstract.

Koala Genome Survey: An Open Data Resource to Improve Conservation Planning.

Genome sequencing is a powerful tool that can inform the management of threatened species. Koalas (Phascolarctos cinereus) are a globally recognized species that captured the hearts and minds of the world during the 2019/2020 Australian megafires. In 2022, koalas were listed as 'Endangered' in Queensland, New South Wales, and the Australian Capital Territory. Populations have declined because of various threats such as land clearing, habitat fragmentation, and disease, all of which are exacerbated by climate change. Here, we present the Koala Genome Survey, an open data resource that was developed after the Australian megafires. A systematic review conducted in 2020 demonstrated that our understanding of genomic diversity within koala populations was scant, with only a handful of SNP studies conducted. Interrogating data showed that only 6 of 49 New South Wales areas of regional koala significance had meaningful genome-wide data, with only 7 locations in Queensland with SNP data and 4 locations in Victoria. In 2021, we launched the Koala Genome Survey to generate resequenced genomes across the Australian east coast. We have publicly released 430 koala genomes (average coverage: 32.25X, range: 11.3-66.8X) on the Amazon Web Services Open Data platform to accelerate research that can inform current and future conservation planning.

Evaluating the genetic consequences of population subdivision as it unfolds and how to best mitigate them: A rare story about koalas.

The genetic consequences of the subdivision of populations are regarded as significant to long-term evolution, and research has shown that the scale and speed at which this is now occurring is critically reducing the adaptive potential of most species which inhabit human-impacted landscapes. Here, we provide a rare and, to our knowledge, the first analysis of this process while it is happening and demonstrate a method of evaluating the effect of mitigation measures such as fauna crossings. We did this by using an extensive genetic data set collected from a koala population which was intensely monitored during the construction of linear transport infrastructure which resulted in the subdivision of their population. First, we found that both allelic richness and effective population size decreased through the process of population subdivision. Second, we predicted the extent to which genetic drift could impact genetic diversity over time and showed that after only 10 generations the resulting two subdivided populations could experience between 12% and 69% loss in genetic diversity. Lastly, using forward simulations we estimated that a minimum of eight koalas would need to disperse from each side of the subdivision per generation to maintain genetic connectivity close to zero but that 16 koalas would ensure that both genetic connectivity and diversity remained unchanged. These results have important consequences for the genetic management of species in human-impacted landscapes by showing which genetic metrics are best to identify immediate loss in genetic diversity and how to evaluate the effectiveness of any mitigation measures.

Human milk oligosaccharides, antimicrobial drugs, and the gut microbiota of term neonates: observations from the KOALA birth cohort study.

The infant gut microbiota affects childhood health. This pioneer microbiota may be vulnerable to antibiotic exposures, but could be supported by prebiotic oligosaccharides found in breast milk and some infant formulas. We sought to characterize the effects of several exposures on the neonatal gut microbiota, including human milk oligosaccharides (HMOs), galacto-oligosaccharides (GOS), and infant/maternal antimicrobial exposures. We profiled the stool microbiota of 1023 one-month-old infants from the KOALA Birth Cohort using 16S rRNA gene amplicon sequencing. We quantified 15 HMOs in breast milk from the mothers of 220 infants, using high-performance liquid chromatography-mass spectrometry. Both breastfeeding and antibiotic exposure decreased gut microbial diversity, but each was associated with contrasting shifts in microbiota composition. Other factors associated with microbiota composition included C-section, homebirth, siblings, and exposure to animals. Neither infant exposure to oral antifungals nor maternal exposure to antibiotics during pregnancy were associated with infant microbiota composition. Four distinct groups of breast milk HMO compositions were evident, corresponding to maternal Secretor status and Lewis group combinations defined by the presence/absence of certain fucosylated HMOs. However, we found the strongest evidence for microbiota associations between two non-fucosylated HMOs: 6'-sialyllactose (6'-SL) and lacto-N-hexaose (LNH), which were associated with lower and higher relative abundances of Bifidobacterium, respectively. Among 111 exclusively formula-fed infants, the GOS-supplemented formula was associated with a lower relative abundance of Clostridium perfringens. In conclusion, the gut microbiota is sensitive to some prebiotic and antibiotic exposures during early infancy and understanding their effects could inform future strategies for safeguarding a health-promoting infant gut microbiota.

Individuality and stability of the koala (Phascolarctos cinereus) faecal microbiota through time.

Gut microbiota studies often rely on a single sample taken per individual, representing a snapshot in time. However, we know that gut microbiota composition in many animals exhibits intra-individual variation over the course of days to months. Such temporal variations can be a confounding factor in studies seeking to compare the gut microbiota of different wild populations, or to assess the impact of medical/veterinary interventions. To date, little is known about the variability of the koala (Phascolarctos cinereus) gut microbiota through time. Here, we characterise the gut microbiota from faecal samples collected at eight timepoints over a month for a captive population of South Australian koalas (n individuals = 7), and monthly over 7 months for a wild population of New South Wales koalas (n individuals = 5). Using 16S rRNA gene sequencing, we found that microbial diversity was stable over the course of days to months. Each koala had a distinct faecal microbiota composition which in the captive koalas was stable across days. The wild koalas showed more variation across months, although each individual still maintained a distinct microbial composition. Per koala, an average of 57 (±16) amplicon sequence variants (ASVs) were detected across all time points; these ASVs accounted for an average of 97% (±1.9%) of the faecal microbial community per koala. The koala faecal microbiota exhibits stability over the course of days to months. Such knowledge will be useful for future studies comparing koala populations and developing microbiota interventions for this regionally endangered marsupial.

Genome Reference Assembly for Bottlenecked Southern Australian Koalas.

Koala populations show marked differences in inbreeding levels and in the presence or absence of the endogenous Koala retrovirus (KoRV). These genetic differences among populations may lead to severe disease impacts threatening koala population viability. In addition, the recent colonization of the koala genome by KoRV provides a unique opportunity to study the process of retroviral adaptation to vertebrate genomes and the impact this has on speciation, genome structure, and function. The genome build described here is from an animal from the bottlenecked Southern population free of endogenous and exogenous KoRV. It provides a more contiguous genome build than the previous koala reference derived from an animal from a more outbred Northern population and is the first koala genome from a KoRV polymerase-free animal.

Susceptibility to a sexually transmitted disease in a wild koala population shows heritable genetic variance but no inbreeding depression.

The koala, one of the most iconic Australian wildlife species, is facing several concomitant threats that are driving population declines. Some threats are well known and have clear methods of prevention (e.g., habitat loss can be reduced with stronger land-clearing control), whereas others are less easily addressed. One of the major current threats to koalas is chlamydial disease, which can have major impacts on individual survival and reproduction rates and can translate into population declines. Effective management strategies for the disease in the wild are currently lacking, and, to date, we know little about the determinants of individual susceptibility to disease. Here, we investigated the genetic basis of variation in susceptibility to chlamydia using one of the most intensively studied wild koala populations. We combined data from veterinary examinations, chlamydia testing, genetic sampling and movement monitoring. Out of our sample of 342 wild koalas, 60 were found to have chlamydia. Using genotype information on 5007 SNPs to investigate the role of genetic variation in determining disease status, we found no evidence of inbreeding depression, but a heritability of 0.11 (95% CI: 0.06-0.23) for the probability that koalas had chlamydia. Heritability of susceptibility to chlamydia could be relevant for future disease management, as it suggests adaptive potential for the population.

Expansion of a retrovirus lineage in the koala genome.

Retroviruses have left their legacy in host genomes over millions of years as endogenous retroviruses (ERVs), and their structure, diversity, and prevalence provide insights into the historical dynamics of retrovirus-host interactions. In bioinformatic analyses of koala (Phascolarctos cinereus) whole-genome sequences, we identify a recently expanded ERV lineage (phaCin-ß) that is related to the New World squirrel monkey retrovirus. This ERV expansion shares many parallels with the ongoing koala retrovirus (KoRV) invasion of the koala genome, including highly similar and mostly intact sequences, and polymorphic ERV loci in the sampled koala population. The recent phaCin-ß ERV colonization of the koala genome appears to predate the current KoRV invasion, but polymorphic ERVs and divergence comparisons between these two lineages predict a currently uncharacterized, possibly still extant, phaCin-ß retrovirus. The genomics approach to ERV-guided discovery of novel retroviruses in host species provides a strong incentive to search for phaCin-ß retroviruses in the Australasian fauna.

A targeted approach to investigating immune genes of an iconic Australian marsupial.

Disease is a contributing factor to the decline of wildlife populations across the globe. Koalas, iconic yet declining Australian marsupials, are predominantly impacted by two pathogens, Chlamydia and koala retrovirus. Chlamydia is an obligate intracellular bacterium and one of the most widespread sexually transmitted infections in humans worldwide. In koalas, Chlamydia infections can present as asymptomatic or can cause a range of ocular and urogenital disease signs, such as conjunctivitis, cystitis and infertility. In this study, we looked at differences in response to Chlamydia in two northern populations of koalas using a targeted gene sequencing of 1209 immune genes in addition to genome-wide reduced representation data. We identified two MHC Class I genes associated with Chlamydia disease progression as well as 25 single nucleotide polymorphisms across 17 genes that were associated with resolution of Chlamydia infection. These genes are involved in the innate immune response (TLR5) and defence (TLR5, IFNγ, SERPINE1, STAT2 and STX4). This study deepens our understanding of the role that genetics plays in disease progression in koalas and leads into future work that will use whole genome resequencing of a larger sample set to investigate in greater detail regions identified in this study. Elucidation of the role of host genetics in disease progression and resolution in koalas will directly contribute to better design of Chlamydia vaccines and management of koala populations which have recently been listed as "endangered."

Subtype distribution and expression of the koala retrovirus in the Japanese zoo koala population.

We investigated the proviral copies and RNA expression in koala retrovirus (KoRV)-infected koalas. To ascertain any variation in viral load by institution, age, sex, or body condition score, we quantified KoRV proviral DNA and RNA loads in captive koalas (n = 37) reared in Japanese zoos. All koalas were positive for KoRV genes (pol, LTRs, and env of KoRV-A) in genomic DNA (gDNA), and 91.89% were positive for the pol gene in RNA. In contrast, the distribution rates of KoRV-B, KoRV-C, KoRV-D, and KoRV-F env genes in gDNA were 94.59%, 27.03%, 67.57%, and 54.05%, respectively. A wide inter-individual variation and/or a significant inter-institutional difference in proviral DNA (p < 0.0001) and RNA (p < 0.001) amounts (copies/103 koala ß-actin copies) were observed in Awaji Farm England Hill Zoo koalas, which were obtained from southern koala populations, suggesting exogenous incorporation of KoRV in these koalas. Significant (p < 0.05) age differences were noted in KoRV RNA load (p < 0.05) and median total RNA load (p < 0.001), with loads higher in younger koalas (joeys and juveniles). Thus, the current study provides the distribution of KoRV subtypes in Japanese zoo koala populations and identifies several additional risk factors (sex, age, and body condition) associated with KoRV expression.

A Possible Link between the Environment and Cryptococcus gattii Nasal Colonisation in Koalas (Phascolarctos cinereus) in the Liverpool Plains, New South Wales.

Cryptococcosis caused by yeasts of the Cryptococcus gattii species complex is an increasingly important mycological disease in humans and other mammals. In Australia, cases of C. gattii-related cryptococcosis are more prevalent in the koala (Phascolarctos cinereus) compared to humans and other animals, likely due to the close association that both C. gattii and koalas have with Eucalyptus species. This provides a cogent opportunity to investigate the epidemiology of spontaneous C. gattii infections in a free-living mammalian host, thereby offering insights into similar infections in humans. This study aimed to establish a link between nasal colonisation by C. gattii in free-ranging koalas and the tree hollows of Eucalyptus species, the key environmental source of the pathogen. We (i) detected and genotyped C. gattii from nine out of 169 free-ranging koalas and representative tree hollows within their home range in the Liverpool Plains, New South Wales, and (ii) examined potential environmental predictors of nasal colonisation in koalas and the presence of C. gattii in tree hollows. Phylogenetic analyses based on multi-locus sequence typing (MLST) revealed that the koalas were most likely colonised by the most abundant C. gattii genotypes found in the Eucalyptus species, or closely related genotypes. Importantly, the likelihood of the presence of C. gattii in tree hollows was correlated with increasing hollow size.

Future-proofing the koala: Synergising genomic and environmental data for effective species management.

Climatic and evolutionary processes are inextricably linked to conservation. Avoiding extinction in rapidly changing environments often depends upon a species' capacity to adapt in the face of extreme selective pressures. Here, we employed exon capture and high-throughput next-generation sequencing to investigate the mechanisms underlying population structure and adaptive genetic variation in the koala (Phascolarctos cinereus), an iconic Australian marsupial that represents a unique conservation challenge because it is not uniformly threatened across its range. An examination of 250 specimens representing 91 wild source locations revealed that five major genetic clusters currently exist on a continental scale. The initial divergence of these clusters appears to have been concordant with the Mid-Brunhes Transition (~430 to 300 kya), a major climatic reorganisation that increased the amplitude of Pleistocene glacial-interglacial cycles. While signatures of polygenic selection and environmental adaptation were detected, strong evidence for repeated, climate-associated range contractions and demographic bottleneck events suggests that geographically isolated refugia may have played a more significant role in the survival of the koala through the Pleistocene glaciation than in situ adaptation. Consequently, the conservation of genome-wide genetic variation must be aligned with the protection of core koala habitat to increase the resilience of vulnerable populations to accelerating anthropogenic threats. Finally, we propose that the five major genetic clusters identified in this study should be accounted for in future koala conservation efforts (e.g., guiding translocations), as existing management divisions in the states of Queensland and New South Wales do not reflect historic or contemporary population structure.

Koala cathelicidin PhciCath5 has antimicrobial activity, including against Chlamydia pecorum.

Devastating fires in Australia over 2019-20 decimated native fauna and flora, including koalas. The resulting population bottleneck, combined with significant loss of habitat, increases the vulnerability of remaining koala populations to threats which include disease. Chlamydia is one disease which causes significant morbidity and mortality in koalas. The predominant pathogenic species, Chlamydia pecorum, causes severe ocular, urogenital and reproductive tract disease. In marsupials, including the koala, gene expansions of an antimicrobial peptide family known as cathelicidins have enabled protection of immunologically naïve pouch young during early development. We propose that koala cathelicidins are active against Chlamydia and other bacteria and fungi. Here we describe ten koala cathelicidins, five of which contained full length coding sequences that were widely expressed in tissues throughout the body. Focusing on these five, we investigate their antimicrobial activity against two koala C. pecorum isolates from distinct serovars; MarsBar and IPTaLE, as well as other bacteria and fungi. One cathelicidin, PhciCath5, inactivated C. pecorum IPTaLE and MarsBar elementary bodies and significantly reduced the number of inclusions compared to the control (p<0.0001). Despite evidence of cathelicidin expression within tissues known to be infected by Chlamydia, natural PhciCath5 concentrations may be inadequate in vivo to prevent or control C. pecorum infections in koalas. PhciCath5 also displayed antimicrobial activity against fungi and Gram negative and positive bacteria, including methicillin-resistant Staphylococcus aureus (MRSA). Electrostatic interactions likely drive PhciCath5 adherence to the pathogen cell membrane, followed by membrane permeabilisation leading to cell death. Activity against E. coli was reduced in the presence of 10% serum and 20% whole blood. Future modification of the PhciCath5 peptide to enhance activity, including in the presence of serum/blood, may provide a novel solution to Chlamydia infection in koalas and other species.

Comparative analysis reveals the species-specific genetic determinants of ACE2 required for SARS-CoV-2 entry.

Coronavirus interaction with its viral receptor is a primary genetic determinant of host range and tissue tropism. SARS-CoV-2 utilizes ACE2 as the receptor to enter host cell in a species-specific manner. We and others have previously shown that ACE2 orthologs from New World monkey, koala and mouse cannot interact with SARS-CoV-2 to mediate viral entry, and this defect can be restored by humanization of the restrictive residues in New World monkey ACE2. To better understand the genetic determinants behind the ability of ACE2 orthologs to support viral entry, we compared koala and mouse ACE2 sequences with that of human and identified the key residues in koala and mouse ACE2 that restrict viral receptor activity. Humanization of these critical residues rendered both koala and mouse ACE2 capable of binding the spike protein and facilitating viral entry. Our study shed more lights into the genetic determinants of ACE2 as the functional receptor of SARS-CoV-2, which facilitates our understanding of viral entry.

Koala methylomes reveal divergent and conserved DNA methylation signatures of X chromosome regulation.

X chromosome inactivation (XCI) mediated by differential DNA methylation between sexes is an iconic example of epigenetic regulation. Although XCI is shared between eutherians and marsupials, the role of DNA methylation in marsupial XCI remains contested. Here, we examine genome-wide signatures of DNA methylation across fives tissues from a male and female koala (Phascolarctos cinereus), and present the first whole-genome, multi-tissue marsupial 'methylome atlas'. Using these novel data, we elucidate divergent versus common features of representative marsupial and eutherian DNA methylation. First, tissue-specific differential DNA methylation in koalas primarily occurs in gene bodies. Second, females show significant global reduction (hypomethylation) of X chromosome DNA methylation compared to males. We show that this pattern is also observed in eutherians. Third, on average, promoter DNA methylation shows little difference between male and female koala X chromosomes, a pattern distinct from that of eutherians. Fourth, the sex-specific DNA methylation landscape upstream of Rsx, the primary lncRNA associated with marsupial XCI, is consistent with the epigenetic regulation of female-specific (and presumably inactive X chromosome-specific) expression. Finally, we use the prominent female X chromosome hypomethylation and classify 98 previously unplaced scaffolds as X-linked, contributing an additional 14.6 Mb (21.5%) to genomic data annotated as the koala X chromosome. Our work demonstrates evolutionarily divergent pathways leading to functionally conserved patterns of XCI in two deep branches of mammals.

Retroviral integrations contribute to elevated host cancer rates during germline invasion.

Repeated retroviral infections of vertebrate germlines have made endogenous retroviruses ubiquitous features of mammalian genomes. However, millions of years of evolution obscure many of the immediate repercussions of retroviral endogenisation on host health. Here we examine retroviral endogenisation during its earliest stages in the koala (Phascolarctos cinereus), a species undergoing germline invasion by koala retrovirus (KoRV) and affected by high cancer prevalence. We characterise KoRV integration sites (IS) in tumour and healthy tissues from 10 koalas, detecting 1002 unique IS, with hotspots of integration occurring in the vicinity of known cancer genes. We find that tumours accumulate novel IS, with proximate genes over-represented for cancer associations. We detect dysregulation of genes containing IS and identify a highly-expressed transduced oncogene. Our data provide insights into the tremendous mutational load suffered by the host during active retroviral germline invasion, a process repeatedly experienced and overcome during the evolution of vertebrate lineages.

High-throughput immunogenetic typing of koalas suggests possible link between MHC alleles and cancers.

Characterizing the allelic diversity within major histocompatibility complex (MHC) genes is an important way of determining the potential genetic resilience of a population to infectious and ecological pressures. For the koala (Phascolarctos cinereus), endemic diseases, anthropogenic factors and climate change are all placing increased pressure on this vulnerable marsupial. To increase the ability of researchers to study MHC genetics in koalas, this study developed and tested a high-throughput immunogenetic profiling methodology for targeting MHC class I UA and UC genes and MHC class II DAB, DBB, DCB and DMB genes in a population of 82 captive koalas. This approach was validated by comparing the determined allelic profiles from 36 koala family units (18 dam-sire-joey units and 18 parent-joey pairs), finding 96% overall congruence within family profiles. Cancers are a significant cause of morbidity in koalas and the risk factors remain undetermined. Our analysis of this captive population revealed several novel MHC alleles, including a potential link between the DBB*03 allele and a risk of developing cancer. This method offers a reliable, high-throughput protocol for expanded study into koala immunogenetics.